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Abstract 

We consider the phenomenological implications of the assumption that 
baryons are systems of three quarks connected through gluon string junction. 
The transfer of baryon number in rapidity space due to the string junction 
propagation is considered in detail. At high energies this process leads to a sig- 
nificant effect on the net baryon production in hN collisions at mid-rapidities. 
The numerical results for midrapidity inclusive densities of different secon- 
daries in the framework of the Quark-Gluon String Model are in reasonable 
agreement with the experimental data. One universal value A ~ 0.25 for the 
strangeness suppression parameter correctly describes the yield ratios of A/p, 
H/A, and The predictions for pp collisions at LHC energies are also 

presented. 



PACS. 25.75.Dw Particle and resonance production 
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1 Introduction 



The Quark-Gluon String Model (QGSM) and the Dual Parton Model (DPM) 
are based on the Dual Topological Unitarization (DTU) and they describe 
quite reasonably many features of high energy production processes in both 
hadron-nucleon and hadron-nucleus collisions [1-6]. High energy interactions 
are considered as taking place via the exchange of one or several Pomerons, 
all elastic and inelastic processes resulting from cutting through or between 
Pomerons [7]. Inclusive spectra of hadrons are related to the corresponding 
fragmentation functions of quarks and diquarks, which are constructed using 
the Reggeon counting rules [8]. 

In the string models, baryons are considered as configurations consisting 
of three connected strings (related to three valence quarks) called string junc- 
tion (SJ) [9-12]. In the processes of secondary production the SJ diffusion 
in rapidity space leads to significant differences in the yields of baryons and 
antibaryons in the midrapidity region even at very high energies [13]. 

A quantitative theoretical description of the baryon number transfer via SJ 
mechanism was suggested in the 90's. The later experimentally observed p/p 
asymmetry at HERA energies was predicted [14] and in [15] it was noted that 
the p/p asymmetry measured at HERA can be obtained by simple extrapola- 
tion of ISR data. 

Important results on the baryon number transfer due to SJ diffusion in 
rapidity space were obtained in [16] and following papers [17-21]. 

In the present paper, we calculate the inclusive densities of different secon- 
daries and compare them with recent RHIC data [22] for pp collisions at \fs = 
200 GeV. The predictions for secondary production at LHC energies are also 
given. 

2 Baryon as 3q + SJ system 

In QCD, the hadrons are composite bound state configurations built up from 
the quark ipi(x),i = 1, ...N c , and gluon G%(x), a — 1, ...,N% — 1, fields. In the 
string models the colour part of a baryon wave function reads as follows (see 
Fig. 1) [9, 11]: 
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Figure 1: Composite structure of a baryon in string models. Quarks are shown by open 
points. 



B = ^ i (x 1 )^ j (x 2 ) i>k{x 3 ) J ljk (x 1 ,x 2 ,x 3 ,x) , (1) 
J ijk {x u x 2 , x 3 , x) = $j,(a:i, x) $^,(x 2 , x) $ k k ,(x 3 , x) ?' j ' k ' , (2) 

i' 

(3) 

i 

where xi, x 2 , x 3 and x are the coordinates of valence quarks and SJ, respectively 
and P(xi,x) represents a path from X\ to x which looks like an open string 
with ends at X\ and x. 

The baryon wave function in Eq.(l) can be defined as a star (or Y) config- 
uration. The Y baryon structure is supported by lattice calculations [23] . 

This picture leads to some general phenomenological predictions. In par- 
ticular, it opens room for exotic states, such as the multiquark bound states, 
4-quark meson and pentaquark [11, 24, 25]. In the case of inclusive reactions 
the baryon number transfer to large rapidity distances in hadrom-nucleon and 
hadron-nucleus reactions can be explained by SJ diffusion. 



®i{x u x) = 



Texp 
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P(xi,x) 
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3 Inclusive spectra of secondary hadrons 
in the Quark— Gluon String Model 



To perform more quantitative predictions a model for multiparticle produc- 
tion has to be adopted. In the present paper we have used the QGSM for the 
numerical calculations. As it was mentioned above, the high energy hadron- 
nucleon collisions are considered in the QGSM as going via the exchange of 
one or several Pomerons. Each Pomeron corresponds to a cylindrical diagram 
(see Fig. 2a), and thus, when cutting a Pomeron, two showers of secondaries 
are produced as it is shown in Fig. 2b. The inclusive spectrum of a secondary 
hadron h is then determined by the convolution of the diquark, valence quark, 
and sea quark distributions u(x,n) in the incident particles with the fragmen- 
tation functions G h (z) of quarks and diquarks into the secondary hadron h. 
Both the diquark and the quark distribution functions depend on the number 
n of cut Pomerons in the considered diagram. 




Figure 2: Cylindrical diagram corresponding to the one-Pomeron exchange contribution 
to elastic pp scattering (a) and the cut of this diagram which determines the contribution 
to the inelastic pp cross section (b). Quarks are shown by solid curves and SJ by dashed 
curves. 

For a nucleon target, the inclusive spectrum of a secondary hadron h has 
the form [1]: 

dn xe da ^2, . 
dy a inel dx F ^ 

where the functions 4>n{ x ) determine the contribution of diagrams with n cut 
Pomerons and w n is the relative weight of this diagram. Here we neglect the 
contribution of diffraction dissociation processes which is very small in the 
midrapidity region. 
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For pp collisions 



4> h J?) = f; q (x + ,n)fZ(x_,n) + fZ(x + ,n)fZ q (x_,n) 

rht„ „\tht 



X± = 2 



\J Arri^/ s + x 2 ± x 



2(n-l)f^x + ,n)f^x.,n) , (5) 

(6) 



where f qq , f q , and f s correspond to the contributions of diquarks, valence 
quarks, and sea quarks, respectively. 

These functions are determined by the convolution of the diquark and quark 
distributions with the fragmentation functions, e.g. 

i 

f q i (x + ,n) = J u g (x 1 ,n)Gg(x + /x 1 )dx 1 . (7) 

x+ 

The diquark and quark distributions, which are normalized to unity, as well 
as the fragmentation functions are determined by Regge intercepts [8]. 

At very high energies both x + and x_ are negligibly small in the midrapidity 
region. In this case all fragmentation functions, which are usually written [8] 
as Gq(z) = a h (l — z) 13 , are constants, 

G h q (x + / Xl ) = a h , (8) 

and lead, in agreement with [26], to 

% = 9H- (s/so)^- 1 - 4 ■ is/so)**™- 1 , (9) 

corresponding to the only one-Pomeron exchange diagram in Fig. 3, which 
is the only diagram contributing to the inclusive density in the central region 
(AGK theorem [7]). The intercept of the supercritical Pomeron a P (0) = 1 + A, 
A = 0.139 [6], is used in the numerical calculations. 

The diagram in Fig. 3 predicts equal inclusive yields for each particle and 
its antiparticle. However, some corrections to the spectra of secondary baryons 
appear for processes which present SJ diffusion in rapidity space. Although 
these corrections would become negligible at energies asymptotically high, they 
result in a significant difference of the baryon and antibaryon yields in the 
midrapidity region for the currently available energy range. Moreover, this 
difference vanishes only very slowly when the energy increases. 
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Figure 3: One-Pomeron-pole diagram determining secondary hadron h production. 

According to [16], we consider three different possibilities to obtain the net 
baryon charge. The first one is the fragmentation of the diquark giving rise to 
a leading baryon (Fig. 4a). A second possibility is to produce a leading meson 
in the first break-up of the string and a baryon in the subsequent break-up 
[8, 27] (Fig. 4b). In these two cases the baryon number transfer is possible only 
for short distances in rapidity. In the third case shown in Fig. 4c both initial 
valence quarks recombine with sea antiquarks into mesons M and a secondary 
baryon is formed by the SJ together with three sea quarks. 

d R d 

c 
c 
c 
c 
c 

a 

Figure 4: QGSM diagrams describing secondary baryon B production by diquark d: initial 
SJ together with two valence quarks and one sea quark (a), initial SJ together with one 
valence quark and two sea quarks (b), and initial SJ together with three sea quarks (c). 

The corresponding fragmentation functions for the secondary baryon B 
production can be written as follows (see [16] for more details): 

G B (z) = a N v qq -z 2 - 5 , (10) 
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G B qs {z) = a N v qs -z 2 (l-z) , (11) 
Gf s (z) = a N ev M -z 1 - asj (l-z) 2 (12) 

for the processes shown in Figs. 4a, 4b, and 4c, respectively and where is 
the normalization parameter, and v qq , v qs , v ss are the relative probabilities for 
different baryons production that can be found by simple quark combinatorics 
[28, 29]. The fraction z of the incident baryon energy carried by the secondary 
baryon decreases from Fig. 4a to Fig. 4c, whereas the mean rapidity gap be- 
tween the incident and secondary baryon increases. The first two processes 
can not contribute to the inclusive spectra in the central region, but the third 
contribution is essential if the value of the intercept of the SJ exchange Regge- 
trajectory, asj, is close to unity. The contribution of the graph in Fig. 4c has 
a coefficient e which determines the small probability of such baryon number 
transfer. 

In [16] the value asj = 0.5 was used. However, for such value of asj 
different values of e were needed for the correct description of the experimental 
data at moderate and high energies. This problem was solved in [17], where 
it was shown with the help of more recent experimental data that all the data 
can be described with the parameter values 

a SJ = 0.9 and e = 0.024. (13) 

It is necessary to note that the process shown in Fig. 4c can be realized very 
naturally in the quark combinatoric approach [28] with the specified probabili- 
ties of a valence quark recombination (fusion) with sea quarks and antiquarks. 

4 Comparison with the experimental data 

The probabilities w n in Eq. (4) are calculated in the frame of Reggeon theory 
[1]. The normalization constants a n (pion production), ax (kaon production), 
ajy (BB pair production), and (baryon production due to SJ diffusion) were 
determined [1, 2, 6] from the experimental data at fixed target energies, where 
the fragmentation functions are not constants. The values of these parameters 
have not been modified for the present calculations, while the values of corre- 
spondent constants for hyperons have been calculated by quark combinatorics 
[28, 29]. For sea quarks we have 

p : n : A + E : 2° : E~ : SI = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 . (14) 
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The ratio S/L determines the strange suppression factor, and 2L + S = 1. 
Usually in soft processes the ratio A = S/L is assumed to be 0.2-0.35. Inside 
this region it should be considered as a free parameter and in the numerical 
calculation we have used the value A = S/L = 0.25 that leads to the best 
agreement with the data [22]. 

The calculated inclusive densities of different secondaries at RHIC, y/s = 
200 GeV, and LHC, y^s = 14 TeV, energies are presented in Table 1, where one 
can see that the agreement of the QGSM calculations with RHIC experimental 
data [22] is reasonably good. 

Table 1 

The QGSM results for midrapidity yields dn/dy (|y| < 0.5) for different 
secondaries at RHIC and LHC energies. The results for e = 0.024 are pre- 
sented only when different from the case e — 0. 



Particle 


RHIC (y/s = 


200 GeV) 


LHC (y/s 


- 14 TeV) 




e = 


e = 0.024 


Experiment [22] 


e = 


e = 0.024 


7T + 


1.27 






2.54 






1.25 






2.54 




K+ 


0.13 




0.14 ±0.01 


0.25 




K~ 


0.12 




0.14 ±0.01 


0.25 




V 


0.0755 


0.0861 




0.177 


0.184 


V 


0.0707 






0.177 




A 


0.0328 


0.0381 


0.0385 ±0.0035 


0.087 


0.0906 


A 


0.0304 




0.0351 ±0.0032 


0.0867 






0.00306 


0.00359 


0.0026 ±0.0009 


0.0108 


0.0112 




0.00298 




0.0029 ±0.001 


0.0108 






0.00020 


0.00025 


* 


0.000902 


0.000934 


n+ 


0.00020 




* 


0.000902 





*dn/dy(Q- + Q+) = 0.00034 ± 0.00019 



The agreement of the QGSM calculations with RHIC experimental data 
[22] is reasonably good. 
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The ratios of p/p production in pp interactions at ^J~s = 200 GeV as the 
functions of rapidity have been calculated in the QGSM with the same param- 
eters used in [20, 21], and they are in reasonable agreement with the experi- 
mental data [30] if the SJ contribution with e = 0.024 is included, while the 
disagreement is evident for the calculation without SJ contribution (i.e. with 
e = 0). It is necessary to note that at asymptotically high energies the ratio 
p/p in the central region is expected to be equal to the unity, so any deviation 
of the p/p ratio from unity has to be explained by some special reason. One 
can see in Table 1 that at the RHIC energies the SJ contribution makes the 
deviation of pp from unity in the midrapidity region about three times bigger 
than in the calculation without SJ contribution. 

The QGSM predicts the deviation of pp ratios from unity due to SJ con- 
tribution on the level of 3-4% accuracy even at the LHC energy. Without SJ 
contribution these ratios are exactly equal to unity. 

The QGSM calculations [17] predict practically equal values of B/B ratios 
in midrapidity region independently on baryon strangeness, what is qualita- 
tively confirmed by the RHIC data on Au-Au collisions [31]. In the case of fl/fl 
production in pp collisions we obtain a non-zero asymmetry (i.e. more Q than 
f2), that is necessary absent in the naive quark model or in all recombination 
models, since both Q and fl have no common valence quarks with the incident 
particles. 

In Fig. 5 we reproduce the experimental data on ratios of yields of different 
secondaries [22] together with our calculations. Agreement is good except 
for only the point of the p/n~ ratio. From the comparison of our results 
with experimental data presented in Table 1 and Fig. 5 we can conclude that 
the universal parameter A = 0.25 describes the ratios of A/p, S/A, and fl/E 
production in a reasonable way. 

5 Conclusion 

We discuss the role of string junction diffusion in the baryon charge transfer 
over large rapidity distances for the cases of pp collisions at RHIC and LHC 
energies. The inclusion of the S J contribution provides a reasonable description 
of the main bulk of the existing experimental data. The calculations of the 
baryon/antibaryon yields and asymmetries without SJ contribution [16, 17] 
clearly diverge for most of the experimental data, where this contribution 
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Figure 5: Ratios of different secondaries produced in midrapidity region in pp collisions at 
sfs — 200 GeV. Short horizontal solid lines show results of the QGSM calculations. 

should be important. Similar results for antibaryon to baryon production 
ratios at RHIC and LHC energies are presented in [33]. They are obtained in 
the framework of DPMJET-III Monte Carlo. Some numerical difference with 
our results comes mainly from the different values of a S j parameter. 
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